Dynamical Breakup of the Fermi Surface in a doped Mott Insulator 
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The evolution from an anomalous metallic phase to a Mott insulator within the two-dimensional 
Hubbard model is investigated by means of the Cellular Dynamical Mean-Field Theory. We show 
that the density-driven Mott metal-insulator transition is approached in a non-uniform way in 
different regions of the momentum space. This gives rise to a breakup of the Fermi surface and to 
the formation of hot and cold regions, whose position depends on the hole or electron like nature of 
the carriers in the system. 
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The mechanism of high-temperature superconductiv- 
ity remains mysterious after nearly two decades of re- 
search. Among many theoretical proposals, a line of 
thought stresses the importance of the proximity to a 
Mott insulating state and of understanding the anoma- 
lous metallic state originating from doping this Mott in- 
sulator, a state which is not described by Fermi-liquid 
theory ; 1J . Along this path two very different approaches 
have been implemented: variational wavefunction studies 
and slave boson approaches . These methods allow 
to connect the Mott insulator to an anomalous normal 
state, interpreting the phenomenon of high temperature 
superconductivity as a direct consequence of the proxim- 
ity to the Mott Metal-Insulator Transition (MIT). Their 
early successes include the prediction of the d-wave sym- 
metry of the order parameter, the existence of a pseudo- 
gap close to half-filling and the dome-like shape of the 
critical temperature. 

In this letter we investigate this issue using the re- 
cently developed Cellular Dynamical Mean-Field Theory 
(CDMFT) |4(|, applied to the two-dimensional Hubbard 
model on the square lattice. The CDMFT method re- 
duces a quantum many body problem to a cluster of sites 
embedded in an effective medium described by a Weiss 
function. This quantity is determined self-consistently 
and can accommodate various broken symmetries, such 
as antiferromagnetism or d-wave superconductivity. In 
this work we consider a metallic phase which does not 
break any symmetry and follow its evolution as a func- 
tion of parameters, such as the doping with respect to 
the Mott insulator. States with long-range order, includ- 
ing superconductivity and antiferromagnetism ||, will 
be discussed in a following paper. CDMFT offers the 
advantage of treating both the coherent and the incoher- 
ent excitations (quasiparticle peak and Hubbard bands) 
on the same footing, capturing the short-range physics 
of singlet formation on links, using self-energies which 
provide a dynamical (i.e. frequency dependent) general- 



ization of the condensates of the early slave boson mean 
field theory. The Hubbard model Hamiltonian is: 
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where Cj CT (c\ a ) are destruction (creation) operators for 
electrons of spin a, is the density of cr-spin elec- 
trons, U is the on-site repulsion and fj, the chemical 
potential which determines the electron density n — 
12icr( c ia c i<?)' is the total number of sites. The 
hopping amplitude ty is chosen to be limited to nearest- 
neighbors t and to next-nearest-neighbors t' , and a strong 
repulsion U = 16t is assumed. In the following we will 
study the hole-doped system (n < 1) with t 1 = — 0.3i, and 
the electron-doped system with t' = +0.3t and n < 1 
(equivalent to the n > 1 system upon a particle-hole 
transformation which reverses the sign of t'). Our main 
results can be summarized as follows: 

i) the normal state approaching the MIT exhibits the 
phenomenon of momentum space differentiation , sim- 
ilar to the results of [7J . Different regions of momentum 
space behave very differently : hot (resp. cold) regions 
are characterized by a small (resp. large) quasiparticle 
residue and a large (resp. small) inverse scattering 
rate ImS. This effect gets most pronounced close to the 
MIT point. 

ii) while the features described in i) are characteristic 
of proximity to the MIT, their detailed mechanism and 
the location in fc-space of the cold and hot regions depend 
on the band parameter t' /t. In particular we find that 
the sign of t'/t for hole-like [electron- like] doped cuprates 
corresponds to a location of the cold region in the fc-space 
around (^, ^) [(0, 7r), (tt, 0)]. This finding agrees with ex- 
perimental observations [6j and earlier theoretical stud- 
ies 0, 0, 0, El 0| • Our study is complementary, but 
consistent, with direct exact diagonalization studies or 
cluster perturbation theory on large clusters Intro- 
ducing the CDMFT self-consistency condition allows for 
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a continuous set of dopings and the high resolution in fre- 
quency needed to study the formation and destruction of 
quasiparticles, as demonstrated in single-site Dynamical 
Mean-Field Theory (DMFT) [l3| and, with remarkable 
success, in tests in one dimension In this paper, 

we argue that the proximity of the MIT is central in this 
phenomenon and we investigate this connection in detail. 

As mentioned above, in CDMFT we have to solve 
for a cluster of sites embedded in a self-consistent non- 
interacting bath. Our cluster is a 2 x 2 plaquette (Fig. H), 
and the bath is truncated to 8 sites, in order to solve the 
impurity model by means of a variant of the Exact Di- 
agonalization (ED) algorithm used in single-site DMFT 
jl5j ] . For the Matsubara frequency calculations we use an 
effective inverse temperature which we set equal to 
/3 = 128 in units of the half bandwidth At. The computa- 
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FIG. 1: Left side: the CDMFT Cluster. Right side: staggered 
magnetization as a function of the next-nearest hopping t' at 
half-filling. The parameter t' controls the magnetic frustra- 
tion in the system. 

tion produces the self-energies E^- on the plaquette and 
the lattice self-energy given by (using the square symme- 
try): 



Elatt(fc, t0 



>) = Eu(w) + Ei 2 (w)(cosfca; -I- cos k y )+ 

E 13 (w) cos k x cos k y (1) 

The cluster self-energies on the imaginary frequency axis 
(where the Green's functions of our ED discrete system 
are smooth jj^j) are presented in Fig. At large dop- 
ing, we recover single-site DMFT results: Ei2,Ei3 ~ 
and only En appreciably different from zero. The lat- 
tice self-energy Ei att (fc,w) JTJ is thereby fc-independent. 
However the cluster self-energies increase sizably at low 
doping as we get close to the MIT and Ei att becomes 
/c-dependent. The zero frequency limit of the real part 
of the self-energy determines the shape of the interact- 
ing Fermi Surface (FS), which is given by t e s(k) = fJL, 
where t g(k) = t{k) — ReEi at t(fc, u — + ), being t(k) the 
Fourier transform of the hopping amplitude Uj . We find 
that the renormalization of the FS becomes appreciable 
close to the MIT. The self-energy itself depends weakly 
on the sign of t', and in particular it has the same sign 
for both positive and negative t' . However, given its large 
magnitude, when combined with t of different signs, it 
produces interacting FS's of very different shape in the 
electron-doped and hole-doped case (see Fig. 0J. This 
can be understood in terms of the renormalized low en- 
ergy hopping coefficients t e g = t — ReEi2(w = + )/2 
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FIG. 2: Real and Imaginary part of the cluster self-energies for 
t'/t = ±0.3. Notice the strong enhancement of the nearest- 
neighbor an next-nearest-neighbor cluster self-energies £12 
and £13 close to the MIT 



and t' eS = t' - ReEi 3 (w = 0+)/4 presented in Fig. 
(where the lu = + limit is extrapolated from the low- 
est Matsubara frequencies). Correlations act to increase 
the value of t e s- ReEi 3 is negative, and it affects t' eS 
quite differently for t'/t = ±0.3 . In fact it increases 
t' eS in absolute value for t' = — 0.3i, in such a way that 
also the ratio (t' eS /t e g)/(t' /t) weakly increases approach- 
ing the MIT, thereby enhancing the hole-like curvature of 
the FS. This has the effect of flattening horizontally the 
Fermi line in the region of the momentum space close to 
(0, 7r) or (77, 0). Because of the enhancement in t e g in the 
vicinity of the Mott transition point, in this region the 
Fermi line is also weakly doping dependent. This behav- 
ior of the antinodal region of cuprates superconductors in 
the normal state has been observed experimentally (lif . 
On the other hand, ReEi 3 decreases t' cS for t' = ±0.3i, 
driving the FS towards a flatter, nested FS as half-filling 
is approached. This is clearly seen in Fig. 0| where the 
spectral function at zero frequency is presented in the 
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FIG. 3: Renormalization of the hopping coefficients and of 
their ratio as a function of density for t'/t — ±0.3. 
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FIG. 4: A(k, uj = + ) in the first quadrant of the Bril- 
louin zone. From the top: in the first row t' — — 0.3t, 
n = 0.73,0.89,0.96, color scale x = 0.28,0.22,0.12; in 
the second t' = ±0.3t, n = 0.70,0.90,0.95, color scale 
x = 0.82,0.34,0.27; in the lowest row t! = ±0.9i, n = 
0.69,0.92,0.96, color scale x = 0.90,0.32,0.22. The white 
dashed line is the FS given by t e a(k) = (i. 



first quadrant of the Brillouin Zone. 

Let us now turn to the formation of hot and cold re- 
gions. For both cases t'/t = ±0.3, ImEn and ImEi3 
have the same negative sign in contrast with the half- 
filled case of earlier studies on an anisotropic model 0. 
While at T = 0, ImEy(w = 0) = 0, the finite effec- 
tive temperature used in the ED calculation gives rise to 



lifetimes which are strongly dependent on the position 
on the FS and reflect the variation of the quasi-particle 
alms,, 



residue Z, 
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of all the imaginary parts of cluster self-energies increase 
as we approach the Mott transition (see Fig. |5J). As they 
sum up differently to form Z// 1 in different regions of the 
/c-space, according to the fc-dependent coefficient multi- 
plying them, Zk may assume different values on different 
points of the FS. This can be already seen in Fig. 0]and 

m 

In the case t' = ±0.3t (electron-doped), as remarked 
above, the FS is flat and crosses (tt/2, 7t/2). There 
ImSi a tt(fc) ~ ImEn. On the other hand, near (tt, 0) 
[(0,7r)], ImEi a tt(fc) — ImEn — ImEi3 . Since ImEii and 
ImEi3 have the same sign, we see that, while the local 
self-energy En tends to decrease the quasiparticle residue 
and increase the inverse scattering rate while approach- 
ing the Mott transition, its effect is counterbalanced by 
the growth of ImEi3, resulting in the formation of a cold 
patch around (0,7r) [(tt, 0)]. For t' = 0.3t instead the 
effect of the real part of the self-energy is to bend the 
Fermi line slightly away from (tt/2, tt/2). The intersec- 
tion of the Fermi line and the zone diagonal now occurs 
at (0.4 tt, 0.4 tt). At that point cosk x ± cosk y ~ 0.6, 
while cosk x cosk y is numerically much smaller (~ 0.1). 
Hence ImE latt ~ ImE u ± 0.6ImEi 2 ± 0.1ImEi 3 . The 
growth of ImEi2 (which has the opposite sign to ImEu 
and ImEi3) reduces the strong correlation effects along 
the diagonal, and produces a cold patch. This "screen- 
ing" of the correlations is larger than the one that takes 
place near the (0, tt) [ (tt, 0) ] region, where the screening 
is given by ImEi3. The modulation in fc-space becomes 
in fact noticeable as soon as ImEi2 > 1.7|ImEi3|. We 
argue that this effect is a precursor to the localization of 
particles that takes place at MIT. The hot quasiparticles 
localize, showing an insulator-like nature, while the cold 
quasiparticles retain their dispersion as the Mott insu- 
lator is approached. This can be seen also in Fig. |SJ 
where we show the spectral weight A(k,uj) (at 4% dop- 
ing for t'/t = —0.3, 5% for t'/t = ±0.3) as a function 
of the energy for different points in the Brillouin zone, 
along the path (0,0) -> {tt,tt) -> (0,tt) -> (0,0). In the 
case t'/t = —0.3 (the upper panel), a quasiparticle peak 
disperses along the (0,0) — ► (tt, tt), where a cold spot is 
formed at the Fermi level; on the contrary the peak is 
less dispersive passing through the region around (0, tt), 
where there is the hot spot. On the other hand, in the 
case t'/t = ±0.3 (lower panel) the quasiparticle peak is 
dispersing more in the vicinity of the region (0,7r) and, 
while dispersing less, evaporates around (tt/2, tt/2). The 
hot-cold spots are in this way switched. 

For t' — +0.3t the Mott insulating state has antiferro- 
magnetic long range order. To investigate the influence 
of antiferromagnetic correlation in determining the evolu- 
tion of the electronic structure described above, we have 
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FIG. 5: Spectral function A(k,u>) as a function of to along the 
path (0,0) — > (7r,7r) — > (0, 7r) — > (7r,7r). In the upper panel 
i'/t = —0.3 and n — 0.96, in the lower panel t'/t — +0.3 and 
n = 0.95. The upper Hubbard band is out of the energy scale 
covered by the figure. 

studied (see Fig. |2J the model with a large t' = +0.9t. 
For this value of t' , in fact, we can rule out any antifer- 
romagnetic long-range order already at half-filling. The 
right side of Fig. ^ displays the staggered magnetization 
of the system at half-filling, m = J2i(~^Y ( n i1 ~ n il)/N c , 
as a function t'/t, from t'/t = -0.9 to t'/t = +0.9: 
a sudden drop of the magnetization is evident around 
\t'/t\ ~ 0.75. The parent insulator state in this case 
is likely to have valence-bond or flux order ^t|- Despite 
the different nature of the parent insulator, the results for 
t' = +Q.9t are similar to those for t' = +0.3t. This, how- 
ever, does not exclude that qualitative differences may 
appear between big and small t' with a better resolution 
in fc-space (i.e. a bigger cluster). 

To summarize, using Cellular DMFT we have studied 
the MIT in the Hubbard model. The approach to the 
insulator is not Fermi-liquid-like, and anisotropic prop- 
erties appear in the momentum space. In particular, 
close to the MIT there is a breakdown of the FS and 
the appearance of cold and hot regions, whose precise lo- 
cation is the result of an interplay of the renormalization 
of the real and the imaginary parts of the self-energy. 
The emergence of these regions is a consequence of the 
proximity of the density-driven Mott transition. In the 
hole-doped case, the quasiparticles survive in the diag- 
onal of the Brillouin zone, near (tt/2, 7r/2). This state, 
which has a fermionic spectrum with pointzeroes, is remi- 
niscent of the flux phases and can be connected smoothly 
to the quasiparticles of the superconducting state. The 



electron-doped case is completely different. On a techni- 
cal level, it is harder to approach the MIT closely, and 
a first-order phase transition may pre-empt a continuous 
approach to the insulating state. Furthermore, the FS 
is renormalized towards nesting, and the quasiparticles 
survive in a small region around (jr, 0) and (0, 7r). These 
quasiparticles cannot be easily deformed into the super- 
conducting state as opposed to the quasiparticles which 
live around (tt/2, tt/2). These properties of the underly- 
ing normal state of the Hubbard model have striking re- 
semblance to what is observed in the cuprates. The hole- 
doped materials have a superconducting region which 
appears almost immediately after doping the Mott in- 
sulator. This superconducting state evolves continuously 
into the pseudogap state which in turn evolves contin- 
uously from the Mott insulator. In the electron-doped 
case, the pseudogap region, if it exists, is small, and a 
much larger doping is needed to reach the superconduct- 
ing phase, which does not seem connected with the Mott 
insulator. 

After the completion of this work, we became aware 
of a related work 0] on the t' = case for the two di- 
mensional Hubbard Model. We thank K.M. Shen and 
Z.-X. Shen for sharing their experimental results, A-M. 
Tremblay, D. Senechal, T. Stanescu, C. Castellani and 
A. Georges for useful discussions. Massimo Capone ac- 
knowledges Italian MIUR Cofin 2003 and the NSF sup- 
port under grant DMR-0096462. 
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